Madelung Energy of the Valence Skipping Compound BaBiC*3 
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Several elements show valence skip fluctuation, for instance, Tl forms the compounds in valence 
states +1 and +3, and Bi forms in +3 and +5 states. This kind of fluctuation gives rise to a 
negative effective attractive interaction and the Kondo effect. In the compounds of valence skipping 
elements, the carrier doping will induce superconductivity with high critical temperature. For 
example, Bai-^K^BiOs shows high T c which is unlikely from the conventional electron-phonon 
mechanism. The reason for the missing of some valence states in such valence skip compounds 
remains a mystery. We have performed the evaluation of the Madelung potential for BaBiC>3, and 
have shown for the first time that charge-ordered state is stabilized if we take into account the 
polarization of the oxygen charge. We argue that the effective Coulomb interaction energy U may 
be negative evaluating the local excitation energy. 

PACS numbers: 74.20.-z, 71.10.Fd, 75.40.Mg 



In the periodic table more than ten elements show va- 
lence skip in all the compounds they form. Tl and In 
form the compounds with the valence states +1 and +3, 
Bi and Sb form with +3 and +5 states, and Pb and Sn 
with +2 and +4 valence states. It is reasonable to expect 
a new microscopic physics behind this phenomenon. The 
reason for missing valence states for these elements still 
remains a mystery[l|,|2,01- The valence skipping leads to 
the Hubbard type model with the negative-?/ attractive 
interaction. Since the attractive interaction possibly pro- 
motes superconductivity, the carrier doping may induce 
superconductivity with relatively high T c . For example, 
Bai-^K^BiOs shows high T c which is unlikely from the 
conventional electron-phonon interaction. A substitu- 
tion of Ba with K for the compound BaBiOj provides 
superconductivity with T c exceeding 30K[1, 01 H 0, 01- 
Bai-^K^BiOs is the perovskite three-dimensional super- 
conductor with clear s-wave gap. The parent compound 
BaBiC>3 is an insulator though BaBiC>3 should be a metal 
according to the band theory 0, [1(3]. A band-calculation 
within the local-density approximation (LDA) leads to 
a (semi)metallic ground state of BaBiO^ [ll[,even under 
the lattice distortion found in the neutron-diffraction ex- 
periments. The average formal valence state of Bi is 
+4. The insulating state of BaBiC>3 would have spin 
1/2 at Bi site if Bi with +4 were allowed. BaBiC>3 
is, however, the insulator with a charge density wave 
(CDW) gap of the order of 2eV. The valence states +3 
and +5 form the CDW state with average charge +4. 
The substitution of K for Ba induces carriers as doped 
holes leading to superconductivity, which is reminiscent 
of the phase diagram of cuprate high temperature super- 
conductors. Since the high critical temperature found 
in Bai-^K^BiOs may be unlikely from the conventional 
electron -phonon mechanism, one can expect that the 
negative-?/ attractive interaction would cooperate in pro- 
moting superconductivity [12, EH G3 • 



The valence skip is important in the study of a charge 
Kondo effect (HI as well as valence- fluctuation induced 
superconductivitv|l6l| . A superconductor Pbi-^TLjTe 
with T c = 1.4K[17| is reasonably expected to belong to 
the same category. The mother compound PbTe is a 
small gap semiconductor. For Tl concentrations up to 
the solubility limit 1.5 percent, the material shows super- 
conductivity with a remarkably high T c as a low carrier 
metal. Tl impurities as negative-C/ centers could leads to 
superconductivity and at the same time Kondo effect as 
a charge analog fla [l8j. 

The purpose of this paper is to investigate the mystery 
of missing valence states evaluating the intra-atomic en- 
ergy for Bi 4+ ion. In general, the total energy is the sum 
of the one-body potential energy and the interaction en- 
ergy. The Coulomb interaction is primarily important 
for the insulator BaBi03. The curvature of the total 
energy as a function of the valence n, which is closely 
related with the intra-atomic energy U, may be presum- 
ably changed due to the Coulomb interaction. In the way 
stated we calculate the Madelung energy of BaBi03. Let 
qi be the charge on the site i, the total Coulomb energy 
is 
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where is the distance between the sites i and j. If wc 



define the potential at the site i as 
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the total Coulomb energy is 
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The potential V$, called the Madelung potential, is eval- 
uated using the Ewald method. (l9l. |2G] We calculated 
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the Madelung potential and the total energy for the va- 
lence skipping compound BaBiC^. We assign the formal 
charges as 



(A) Ba : 2+, Bi : 4+, O : 2- 



(4) 



(B) Ba : 2+, Bi(l) : 3+, Bi(2) : 5+, O : 2- (5) 

where in the latter case the charge density wave with va- 
lences 3+ and 5+ is assumed. Let us denote the number 
of Ba atoms as Na, then the number of Bi and O atoms 
are Na and 3Na, respectively. The total Madelung en- 
ergy for the case (A) is 

E M adei = ^(N A q Ba V Ba + N A q B iV Bi + 3N A qoV ) 



N A 



(qBaV B a + qBiV B i + MoVq), (6) 



where q Ba = 2, q B i = 4 and qo = —2. For the case (B) 
the total Madelung energy is 

EAIadel = (?Ba^Ba + « QBi(l) Vbi(1) + ~ QBi(2) V B i(2) 
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for q Ba = 2, q B i(i) = 3, q Bl ( 2 ) = 5 and q Q = -2. The 
intra-atomic repulsion energy U can be calculated if the 
total energy is obtained. For the valence n, the energy U 
is defined as 



U n — E; 



n+1 



E n -i — 2E n 



(8) 



where Ei is the energy of the valence i state. The U for 
Bi + per two bismuth atoms has been estimated from the 
measured ionization energies of the elements aspj 



U i+ = E ion {Bi^)+E ion {Be + )-2E ion {Be + ) = 10.7eV. 

(9) 

Here Ei on (Bi n+ ) is the ionization energy of the clement 
Bi n+ . This estimated value is quite large although it is 
smaller than for Bi 3 + (f7 3+ = 19.7eV) and Bi 5 + (U 5+ = 
32.3eV). Thus the simple estimate results in the positive 
U. 

We compare the total energy of the uniform state (A) 
and the ordered CDW state (B), where the total energy is 
the sum of the total Madelung energy and the one-body 
(ionic) potential energy . The total energy difference is 

AE = E{Ba\ + Bi 3+ Bi 5+ Ol-)- E{Ba 2 +Bi i+ Bi i+ Ol-) 
+ E lon (Bi 5+ )+E lon (Bi 3+ )-2E lon (Bi A+ ), (10) 



where 



E(Bal + Bi 3+ Bi 5+ Ol~) 



and 



E{Ba\ + Bi 4+ Bi i+ Ol~) are the total Madelung energies 
per two BaBiOs's, respectively. The charge-ordered 
state is stabilized if AE is negative. The Madelung 
potential V a (a=Ba,Bi and O) are shown in the Table 
I where we use the convention that the potential Vt 
is negative for the positive charge qi > 0. The total 
potential energy per BaBiC>3, E tot /NA, is -81.95 eV for 



(A) and -83.40 eV for (B), respectively, which are shown 
in the first and second rows of the Table I. Apparently 
the case (B) with the charge density wave has the lower 
energy. From these values, we obtain AE — 7.8eV 
which is still positive. It is now important to notice that 
the AE is reduced if we consider the movement of the 
oxygen ions. The neutron diffraction demonstrated that 
the Bi ions occupy two crystallographically inequivalent 
sites with different Bi-0 bond lengths[2l[. It has been 
reported by the neutron diffraction experiment that 
the oxygen atoms are out of the center of two Bi ions 
by 4 percent of the Bi-0 bond length due to lattice 
relaxation [io| (see Fig.l). Thus it is reasonable to move 
the positions of oxygen atoms with different Bi-0 bond 
lengths in calculating the Madelung potential. In order 
to take into account the polarization of electrons of 
oxygen atoms, we have further moved the center of point 
charges more than reported by the neutron diffraction. 
It is plausible to assume that the center of the negative 
charges is more close to the neighboring positive charges 
Bi 5+ than the nucleus of the atom. The evaluated 
results are shown in the Table I, and are also shown 
in Fig. 2 as a function of the bond length. We found 
that the stabilization energy AE is reduced linearly 
as a function of the bond length O-Bi(l) and becomes 
negative for about 10 percent movement. This indeed 
shows the stability of the observed CDW state. 

Moreover, chemical bonding between Bi and O ions 
can further stabilize this CDW state. Bi 5+ ion can gain 
energy as the neighboring oxygen ions move toward it, 
lowering the energy of the oxygen states while raising 
the energy of the empty Bi s states, as compared to the 
Bi 3+ where those s-states are occupied. We consider Bi 5+ 
Og~ cluster and the three-parameter model shown by 
Mattheiss and Hamann[ll| for simplicity. Adopting more 
precise five-parameter model does not change the result 
seriously. The relevant orbitals are \Bi — s > and \0 — 
Alg >= (\pi > +\P2 > +\P3 > +\Pi > +\P5 > +|P6 > 
)/V6 orbitals, where \p n > (n — 1, ..,6) denotes the n- 
th oxygen p-orbital directing the Bi site. The secular 
equation is: 



A - e t 
t -e 



(11) 



here A is the difference of the one-electron orbital energy 
A = e(Bi — s) — e(0 — Ai g ), and t is the effective trans- 
fer matrix t = V&(spa). Using the well-known relation 
(spa) oc d~ 2 [2|| and the distortion x = 0.26, then the 
energy gain (loss) by the distortion is 0.42eV per electron 
for the bonding (antibonding) state. Since the Bi 5+ site 
has two empty states, thus the total energy gain is 0.84eV 
per Bi 5+ .On the other hand, Bi 3+ ion cannot gain energy 
because there is no empty s-states. This effect apparently 
promotes forming the CDW state. 

In the above discussion we only considered the insulat- 
ing states. However, in the uniform BaBiOa the system 
may be metallic. In fact, band-calculations show that 
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uniform BaBiOa with no distortion and oxygen polariza- 
tion should be metallic [Til. [2^ |. The width of Bi6s band 
is about 4eV, and the Fermi level is just the middle of 
this band. Thus considering the small density of states at 
the bottom of this band, the average energy gain forming 
Bloch state is less than leV per electron. Each Bi 4+ has 
one Bi6s electron, so that the total energy of the uniform 
state may go down less than 2eV. However, this value 
is not so large compared with the large energy gain of 
CDW state, and a further but little change of the neg- 
ative charge, namely x ~ 0.285, can compensate this 
effect. 

In order to further investigate the possibility of the 
negative U originating from the long-range Coulomb in- 
teraction, we consider a state excited locally from the 
CDW ground state: a (3+, 5+) Bi pair (e.g. Bi atoms 
at (0,0) and (0.5,0) in Fig.l) is changed into a (4+, 4+) 
pair in the CDW background. The excitation energy for 
this process is 

E loc = \{V B a - V Bi2 ) ft( t ~ i * , i ) ~ Wfi 

Here = a%e x denotes the position vector of Bi(2) at 
e x = (0.5,0,0), and a\ is twice the average Bi-Bi dis- 
tance (ai = 8.700A). Vbh and Vsa are the Madelung 
potential at the Bi(l) and Bi(2) site, respectively. The 
first term in eq.(|12|) is the change of the Madelung en- 
ergy of the Bi pair, and the second term denotes that 
of the Madelung potential of the other ions. The third 
term shows the energy gain of the ionic (one-body poten- 
tial) energy for (3+, 5+) — * (4+, 4+), shown in eq.(9). 
We performed the three-dimensional summation in real 
space for sites in the sphere of the radius R using the 
procedure of Harrison [22|. This summation converged 
well for R ~ 200ao (see Ref.22). We do not suffer from 
the usual extreme oscillation of the potential in this case 
because we only consider the difference of the Madelung 
energy. Thus we do not need to add the —Q/R term in 
contrast to the method in Ref.22. In the case of the Bi(l)- 
O bond-length x = 0.282, the local excitation energy is 
E loc = 6.92 + 10.56 - 10.7 = 6.78eV. Hence the CDW 
state is stable against local valence fluctuation. Rice and 
Sneddon examined the low-energy effective Hamiltonian 
for the CDW state of BaBi03 and insisted that there are 
two types of excitations [23]. One is the single-particle 
excitation mentioned above, and the other is the two- 
particle excitation which is an exchange of local charges 
as (3+, 5+)^ (5+, 3+). The excitation energy for (3+, 
5+) — > (5+, 3+) is extremely large because the second 



term of eq.(ll) is twice as large as the process (3+, 5+) 
— > (4+, 4+). From the definition of U n in eq.©, the 
quantity XJ l £? al = —Ei oca i resembles the intra-atomic en- 
ergy. The value U^f al = — 6.78eV seems too large as 
the 'negative' U ; the high dielectric constant and doped 
carriers, however, may reduce Ei oc . 

In summary, we have evaluated the Madelung poten- 
tial for the valence skipping compound BaBi03 using the 
Ewald method. The total energy of the charge-ordered 

TABLE I: Madelung potential in units of eV. The total elec- 
trostatic energy Em ad/ N a and the stabilization energy of the 
CDW state A_E defined as eq.(10) are shown in units of eV. 
Note that EMad/NA is for one BaBiC^. The bond length be- 
tween bismuth (1) and oxygen atoms is normalized as 0.25 
when there is no lattice relaxation. 



Bond lenj 


;th 


Ba 


Bi(l) 


Bi(2) 


O 


E M ad/N A 


AE 


0.25 


(A) 


-8.917 


-20.487 


-20.487 


10.685 


-81.95 




0.25 


(B) 


-8.917 


-17.594 


-23.379 


10.685 


-83.39 


7.8 


0.26 


(A) 


-8.908 


-19.321 


-21.789 


10.661 


-82.00 




0.26 


(B) 


-8.908 


-16.428 


-24.681 


10.867 


-84.68 


5.3 


0.27 


(A) 


-8.881 


-18.277 


-23.244 


10.589 


-82.17 




0.27 


(B) 


-8.881 


-15.384 


-26.136 


11.003 


-86.10 


2.8 


0.28 


(A) 


-8.838 


-17.341 


-24.874 


10.465 


-82.45 




0.28 


(B) 


-8.838 


-14.448 


-27.766 


11.093 


-87.66 


0.26 


0.282 


(A) 


-8.827 


-17.166 


-25.223 


10.434 


-82.52 




0.282 


(B) 


-8.827 


-14.273 


-28.115 


11.105 


-87.99 


-0.24 



CDW state was calculated using the Madelung poten- 
tial. We have evaluated the excitation energy for the 
local disorder (3+, 5+) — > (4+, 4+). It was shown that 
the polarization of electrons reduces the energy difference 
AE and that an attractive 'negative' U is presumably 
realized if the center of the negative charges is moved 
from the center of the Bi-Bi bond by about 10 percent 
of the bond length. This indicates a possibility that 
the high critical temperature of Bai-^K^BiOa is due to 
the valence-skipping induced negative U as well as the 
electron-phonon interaction. The results also suggests 
that the elements with high electronic polarizability have 
strong possibility to induce negative U. Since the Te 2_ 
ion has high polarizability a negative U induced su- 
perconductivity is likely present in Tl doped PbTe. The 
valence skip may provide a new idea for the material de- 
sign of new superconductors. 
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FIG. 1: The positions of bismuth and oxygen atoms in the 
plane of z = 0. Here we ignored the small tilting of BiOe 
octahedra. 



[9] A. W. Sleight, J. J. Gillson, and P. E. Bierstedt: Solid 
State Commun. 17 (1975) 27. 
[10] C. Chaillout, J. P. Remeika and A. Santoro: Solid State 



Commun. 65 (1988) 1363. 
[11] L. F. Mattheiss and D. R. Hamann: Phys. Rev. B28 
(1983) 4227. 

[12] C. S. Ting, D. N. Talwar, and K. L. Ngai: Phys. Rev. 

Lett. 45 (1980) 1213. 
[13] H.-B. Schuttler, M. Jarrell, and D. J. Scalapino: Phys. 

Rev. B39 (1989) 6501. 
[14] A. G. Mal'shukov: Solid State Commun. 77 (1991) 57. 
[15] A. Taraphder and P. Coleman: Phys. Rev. Lett. 66 

(1991) 2814. 

[16] Y. Ohnishi and K. Miyake; J. Phys. Soc. Jpn. 69 (2000) 
3955. 

[17] Y. Matsushita, H. Bluhm, T. H. Geballe, and I. R. Fisher: 

Phys. Rev. Lett. 94 (2005) 157002. 
[18] M. Dzero and J. Schmalian: Phys. Rev. Lett. 94 (2005) 

157003. 

[19] C. Kittel: Introduction to Solid State Physics 8th ed. 
(New Jersey, John Wiley and Sons Inc., 2005). pp. 644- 
648. 

[20] J. Kondo: J. Phys. Soc. Jpn. 58 (1989) 2884. 
[21] E. E. Cox and A. W. Sleight: Solid State Commun. 90 
(1976) 969. 

[22] W. A. Harrison: Phys. Rev. B73 (2006) 212103. 
[23] T. M. Rice and L. Sneddon: Phys. Rev. Lett. 47 (1981) 
689. 

[24] K. Takegahara and T. Kasuya: J. Phys. Soc. Jpn. 56 
(1987) 1478. 

[25] W. A. Harrison: Elementary Electronic Structure (World 
Scientific, Singapore, 1999), revised edition (2004). 





-164F 
-166 - 
-168 - 
-170 - 
-172 - 
-174 - 
-176L 



FIG. 2: The Madelung potential at the Ba, Bi and O sites 
and the total energy as a function of the bond length. The 
solid lines and the dotted lines correspond to the case (A) 
and the case (B) (with the charge density wave), respectively. 
The bond length is normalized as 0.25 when there is no lattice 
relaxation. 



